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ABSTRACT 



A study was conducted of the amplitude modulation of an 
acoustic signal forward scattered from a wind-generated, 
model ocean surface. The demodulated spectrum of the surface 
scattered sound, was compared with the model ocean spectrum 
while varying the acoustic frequency and angle of incidence. 
Several methods of enhancing signal detection through a 
knowledge of the ocean spectrum were studied. The most 
promising of these shows potential gains of the order of 
three decibels when the received signal is weak and 
fluctuating greatly, as at a Lloyd Mirror minimum. 
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I. INTRODUCTION 



It has long been known that received signals from 
distant underwater sound sources vary In time. The effect 
has been described as "builds and fades" of the signal. 

Though this effect is automatically interpreted by the 
human ear, which listens when the sound is present and tunes 
out when it disappears in the noise, little or no advantage 
of the phenomenon has been taken by the extensive signal 
processing equipment available today. Though nominally far 
superior to the ear in response-time and frequency resolution, 
this equipment often fails to achieve even parity with a 
human observer when processing a time varying, or modulated, 
signal. This is due primarily to the electronic equipment's 
narrow bandwidth and the fact that it integrates continuously 
rather than selectively. 

It was the objective of this research to investigate some 
of the parameters affecting this phenomenon of a modulated 
acoustic signal with particular attention to determining, for 
surface affected signals: 

1. The best laboratory modeling techniques. 

2. The relation of the modulation frequency to the 
surface wave frequency. 

3. The variation of the modulation amplitude and spectrum 
with frequency and source/receiver directivity. 

The optimization of signal processing of an amplitude 
modulated signal by digital techniques. 
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In short an attempt was made to discover a ’’temporal 
filtering" scheme which could be used in concert with 
existing high speed frequency filtering techniques, such as 
the Past Fourier transform, to enhance signal processing. 

The concept has been dubbed "EASE", Environmentally Adjusted 
Signal Enhancement. 



II. THEORY 



The subject of scattering of sound from the ocean 
surface has been treated by several authors utilizing nearly 
as many methods. The following is a summary of Refs.l and 2 
which point up those features of the scattered spectrum 
which are germane to understanding the relationships between 
sound fluctuations in the upper ocean and the ocean wave 
spectrum. 

As Parkins summarizes the phenomenon: "When a mono- 

chromatic acoustic wave is incident on the surface of the 
ocean the envelope of reradiation varies in time or fades." 

Assuming that the ocean is generated by a stationary 
homogeneous, Gaussian process, and that it is smooth enough 
so that multiple scattering can be neglected, the incident 
acoustic wave of harmonic time variation e^ 10 ^ is given by: 



*i - 



.+ - S l' ? lwt 
$ e e 
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k l = k lx* + k ly y + k lz Z 



2 TT 



= — [2.x + my + nz] 



4> + is an amplitude term, l, m, and n are the direction 
cosines and X is the wave length of the incident acoustic 
wave. Other factors are defined In Figure 1. 
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FIGURE 1. Coordinate System for Theoretical Development 

Following the development of Beckmann and Spizzichino 
[Ref. 2] which stipulates a source free, isovelocity ocean, 
and using the Helmholtz Integral development of Stratton 
[Ref. 3]> the farfield reradiated pressure may be expressed 
as : 
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with; 



^ = K x x + K y y + K z z = k Q -k x 

and other factors again as defined in Figure 1. 

The irradiated area has been assumed to be an LxL 
aperture which is several wavelengths across. 

In this research it was the spectral density of the 
reradiated pressure that was of interest. This may be 
obtained by first forming the autocovariance function, 

r(t) = <$ r (x ,y ,t ) $ r *(x,y,t+T)> 

The autocovariance function having been found, the Wiener- 
Khinchin Theorem states : 

S(ft) =5 T{R(t)} = I R(t) e lfiT dt 

_oo 

where ft is the frequency of the reradiated (or scattered) 
sound. 

Continuing with Parkins, and assuming that the surface 
irregularities are much smaller than the aperture, the 
expression for the spectral density of the scattered sound 
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is as given below. It should be noted at this point that 
the requirement for small surface irregularities corresponds 
to a roughness parameter. 



•yfg = K z a = [cos 0-^ + cos 0 2 ] 

as defined by Beckmann and Spizzichino, much less than one. 
a is the ocean RMS wave height , X the acoustic wave length 
and 0-^ and © 2 the angles of incidence and reflection, 
respectively (see Section IV. A. 4 for a further discussion 

of-^/g"). ' 

The expression for the spectral density of the scattered 
sound is: 

-K 2 a 2 16 sin 2 K Z sin 2 K m 

S(fi) = Be z [6(n-u) 5— — 5 £_ 



+ 2£ 2 K 2 a 2 { [A* (K ,K )] 2 6[fi-co-(K 2 +K 2 ) 1< (G) Js ] 

L A j A y 

+ [A'(-K y ,-K y ) ] 2 6[fi-a)+(K x 2 +K y 2 ) !s (G) J2 ]}] 



with: G the acceleration of gravity 

A* the Neumann-Pierson directional wave spectrum 
B the amplitude term 
6 the Dirac delta 
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Parkins summarizes this equation: "....(this equation) 

shows the signal reradiated by a wind generated sea to have 
a specularly reflected component occurring at the incident 
frequency w and two scattered components which are Doppler 
shifted by equal amounts from co." 

The spectrum of the reradiated sound is thus seen to be 
frequency modulated by the Doppler shift term and amplitude 
modulated by the wave spectrum. It is the amplitude 
modulation by the ocean waves that is considered in this 
research. 
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Ill . EXPERIMENTAL SETUP 

■i n ■ ■ ■ *— - tf * 1 

This research was carried out using the Naval Postgraduate 
School OCEAN ACOUSTIC WAVE FACILITY (OAWF) with data 
collection and analysis accomplished with a hybrid computer 
setup designed by the Special Projects Group of NADC and 
affectionately referred to as OPHELEA (OCEAN PHYSICS 
ENVIRONMENT EFFECTS ANALYZER). All other equipment utilized 
was standard commercial electronic gear and the author feels 
that a detailed explanation of its functions and capabilities 
is unnecessary (see Equipment List - Section III.D). 

A . OAWF 

This facility consists of the unique combination of a 
wind-water tunnel and an anechoic tank (see Figures 2 and 3)* 
Up to five 3/4 HP centrifugal fans (three in this case) 
comprise a Sea Exciter (Figure 4) which forces air down the 
wind-water tunnel, which is about 17 meters in length, 
generating waves which empty into an anechoic tank. The air 
space in the tunnel was standardized at 14.25 cm and the 
cross section of the entire tunnel measures about 1.2 meters 
square. Observation of the waves at various points along 
the tunnel is permitted by glass observation windows shown 
in Figure 5. Waves entering the 3x3x2. 9 meter tank portion 
of the facility are constrained from diffracting by a 
plexiglass restraining barrier and are terminated by a 
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"beach" to prevent reflections. The beach spans the width 
of the tank experimental area and consists of a screen 
covered wedge filled with stainless steel millings which 
cause incident waves to dissipate rather than reflect. 

The tank itself is lined with redwood 4x4* s on all sides 
and the bottom which ensures excellent sound absorption at 
the boundaries as pointed out in Ref. 4. These redwood 
beams are oriented with a corner facing toward the tank 
interior, except on the wall bordering the waves which is 
flat, to further increase absorption. A movable race permits 
two-dimensional positioning of transducers to within one 
centimeter. This facility models the real ocean on a scale 
of about 50-1 and is shown in operation in Figures 6 through 
8 . 



B. OPHELEA 

The OPHELEA system consists of essentially three 
components which are interfaced to provide high speed analog 
to digital data collection and digital data processing with 
somewhat slower display. The design was developed by the 
Special Projects Section of the Naval Air Development Center 
in conjunction with Pinkerton Computer Consultants Inc. , 
of Warminster, Pennsylvania. The three components are shown 
in Figure 9 and described below. 

1. Interdata Model 70 Computer 

This minicomputer is modularly constructed and 
features LSI and other state-of-the-art computer technology. 
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Figure 6. OAWF, Looking Down into Anechoic Tank 
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Figure 7. OAWF in Operation with Wave Peaks Shown Against the 
Plexiglass Restraining Barrier 
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Figure 8. OAWF, Looking from the Tank Back up the Wind/Wave Tunnel 
Toward the Sea Exciter 




Figure 9. The OPHELEA System with Computer and A/D Converter 
in Large Cabinet to the Left 
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It is a 1 6 Bit half-word oriented digital computer that is 
FORTRAN programmable and has sufficient memory (32 thousand 
8-bit bytes) capacity to process in excess of 8000 data 
points with an average size program. In this research it 
was programmed to perform a Fast Fourier Transform on input 
data and output spectral density vs. frequency. Some 
statistical analysis was also performed. This minicomputer 
occupies the upper half of the large cabinet shown in 
Figure 9 . 

2 . Phoenix Analog to Digital Converters Model ADC 712 
Two ADC 712 converters may be used separately or 
simultaneously. Each A/D converter is a high speed device 
capable of encoding up to ±10 volt input signals into 
12 bits (11 bits + sign) of data, providing resolution to 
one part in ^095 (accuracy of 5mv in 20v). The maximum rate 
of conversion is 6.5 microseconds and it has performed 
reliably at a sampling rate of up to 300kHz. The A/D is 
driven by a "command to convert" signal consisting of a 
5 volt positive pulse of at least 200 microseconds duration 
which is supplied by an external oscillator. Measurement 
is accomplished by successive comparisons against an internal 
reference voltage after which the value is converted to a 
digital number. The A/D converters are located in part of 
the lower half of the OPHELEA cabinet. 
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3. Texas Instruments Silent 700 Electronic 

Data Terminal Model 733 

The TI 733 consists of a keyboard which is used as 
a programming and input/output control device, a printer 
which outputs data or answers queries of the computer on 
heat sensitive paper and a play back and record section 
which is used to read in computer programs and record/ 
playback data. In this research, due to storage limitations, 
data were sometimes taken A/D, read out to tape cassettes, 
then reread in for processing after the appropriate program 
had been loaded. A capability of interfacing with the 
IBM- 360/67 CP/CMS system, allowing direct access to the NPS 
computer facility, was also utilized in processing some 
statistical wave height information. 

C. SEQUENCE OP OPERATION 

The sequence of operation used most often in this 
research was as follows: 

1. An analog signal (varying sound pressure, instan- 
taneous wave height, etc.) was input to the A/D converter 
which, being triggered at twice the Nyquist rate (a value 
found to result in minimal aliasing effects), converted it 
to digital data. 

2. Those data were then either: 

a. directly processed by the appropriate program 
which had been installed in the computer, or, 

b. read out to tape cassettes for later analysis 
as described above. 
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3. Reduced data and desired results were displayed on 
the paper readout of the terminal. Figure 10 is a pictorial 
description of this process. 

D. EQUIPMENT LIST 

In order to avoid needless repetition of the lengthy 
manufacturers’ names for electronic equipment throughout 
this work, a number of abbreviations will be used instead. 
These short titles will be employed for the most part in 
block diagrams and schematics but may also appear in the 
text. The abbreviations, with the corresponding full 
equipment title/description are as listed below: 



ABBREVIATION 
CORRELATION COMPUTER 

DC POWER SUPPLY 

G.R. OSCILLATOR 

HP PRECISION AMP 

HP POWER AMP 

KH FILTER- 3.3 ^2 
KH FILTER- 35 50 
O’ SCOPE 



MANUFACTURERS TITLE 



Princeton Applied Research 
Model 101 Correlation Function 
Computer 

Hewlett-Packard Model 721A DC 
Power Supply 

General Radio Model 1312 Decade 
Oscillator 

Hewlett-Packard Model 463A 
Precision Amplifier 

Hewlett-Packard Model 
Power Amplifier 

Khron-Hite Model 33^2 Filter 

Khron-Hite Model 3350 Filter 

Tektronix, Inc. Model 5^5B 
Oscilloscope 
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ABBREVIATION 


MANUFACTURERS TITLE 


PAR PREAMP 


Princeton Applied Research Model 113 
PreAmplifier (Note: HFRO and LFRO 

refer to the High/Low frequency roll 
off feature of the preamp — a 
filtering option) 


WAVETEK 142 


Wavetek Model 142 HF VCG Generator 


WAVETEK 144 


Wavetek Model 144 HF VCG Generator 


X-Y RECORDER 


Varian Aerograph Model F-100 Strip 
Chart Recorder 
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DATA 

(digital 

samples) 



DATA 



PRINTED 

OUTPUT 



Figure 10. Data Processing Sequence 
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IV. EXPERIMENTAL PROCEDURES AND RESULTS 
A. OAWF CHARACTERISTICS 

In order both to ensure the validity of the research to 
follow and to gain an insight into the significant parameters 
affecting the scattering of sound from a real or model ocean, 
several preliminary experiments were conducted. In the 
sections that follow each will be fully discussed with setup 
diagrams, procedural explanations and results included. 

1 . Frequency Spectrum 

In determining the frequency spectrum of the ocean 
wave heights, a highly linear probe designed at the Civil 
Engineering Department of Stanford University and pictured 
in Figures 11 and 12 was utilized. The probe is capacitive 
in nature, with a covered wire as the sensor. The wire 
serves as one electrode, the water as the other and the 
insulation is the dielectric of the capacitive element. 

The probe is excited by a 2.4 kHz sinusoidal carrier 
(amplitude 4.5 volts peak to peak). When the probe is 
immersed in the water wave system the output across a 
capacitive bridge (Figure 13) is amplitude modulated by 
the passing waves rising and falling around the wire and 
thereby changing the effective length of the capacitor. 

This AM signal is demodulated by the detector shown in 
Figure 1 4 and the resultant time varying DC signal, a replica 
of the instantaneous wave height at one point, is frequency 
analyzed by the OPHELEA system. 
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Figure 12. 

Wave Probe in 
Operation 
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Figure 13. Wave Height Capacitive Bridge 



TA-38 IN271 





Figure 1 J I. Wave Height Detector (Demodulator) 
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As noted above frequency analysis was accomplished 
by the OPHELEA system, using a program called PING 1C, a 
copy of which is included in Appendix A. In this case the 
program called for taking 256 A/D samples at a sampling rate 
of 20 Hz and computing a 128-point spectral density versus 
frequency output. These parameters provide a frequency 
resolution of greater than .08 Hz. 

When taking A/D samples at sampling frequencies in 
excess of 20 kHz, a sinusoidal oscillator which is half-wave 
rectified by computer circuitry was used to provide the 
"command-to-convert " pulse train. At low frequencies however 
the sinusoidal pulse duration is too great and erroneous 
multiple triggering of the A/D converter results. Therefore 
the WAVETEK 1*14, with a positive pulse waveform and with its 
asymmetry control set for minimum pulse duration, was employed 
for the low sampling frequencies required for finding spectra 
of wave and demodulated scattered sound. A block diagram 
of the experimental setup is shown in Figure 15 and the 
normalized spectral density plot is shown in Figure 16. 

This plot is the average of forty runs conducted as outlined 
above (i.e., *10-128 point FFT analyses). 

As can be seen from the figure, the predominant OAWF 
ocean frequency is about 2.5Hz as compared with a typical 
actual ocean frequency on the order of 0.1Hz. Though neither 
ocean is monochromatic, the OAWF spectrum is fairly narrow. 

Its bandwidth is about 0.6Hz at the half power points and 
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GAIN = 5K 
HFRO = 30 Hz 
LFRO = DC 




10 Hz 



WAVETEK 

144 



20 Hz 




TO OPHELEA 
■^(Sampling Freq.) 



* See Figure 13 

# See Figure 14 

Figure 15. Experimental Setup for Determining OAV/F 
Surface Wave Spectrum 
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Figure 16 . Normalized Wave Height Frequency Spectrum (Average of 40 Runs) 



the spectral density tends to fall off less rapidly toward 
higher frequencies than does the actual ocean spectrum 
[Ref. 5]. This is due to the greater contribution of 
capillary waves in the OAWF waves. The irregularity at 
about 2. 4Hz was felt at first to be merely a coincidental 
artifact which would disappear with additional averaging. 
However, though it smoothed some, it has persisted. The 
source of the irregularity may have been lack of adequate 
fetch for w ave growth to an equilibrium spectrum. It is 
noted that a similar effect was observed in the work of 
Fowler and Scheibel [Ref. 6] at short fetch. 

2 . Probability Density Function 

As Kinsman has noted [Ref. 7] the probability density 
function associated with wind driven ocean wave heights is 
nearly Gaussian. In order to verify that the OAWF sea met 
this criterion, wave height samples were collected using the 
setup shown in Figure 15 of the previous section and analysis 
was accomplished using the standard IBM- 360/67 Library 
program "HISTF" . This process was accomplished as follows: 

a. Eight groups of 1024 wave height samples were 
taken from the previously described wave probe using the 
OPHELEA system and a program modestly entitled PERK 1A and 
were stored in the minicomputer memory. 

b. The data were then read out to tape cassettes in 
a revised format compatible with the IBM-360 CP/CMS system. 

c. Data were converted to punch cards using the 
G’P/CMS system and analyzed using "HISTF". 



38 



The "HISTF" program outputs various statistical 
parameters and plots a histogram of the input data. Each of 
the eight data groups was analyzed and the results averaged. 

A graph of the resulting average histogram, as well as a 
comparable normalized Gaussian PDF, is shown in Figure 17. 

As can be seen, the histogram distribution is slightly skewed 
and displays negative kurtosis. Though not exactly agreeing 
with any one of Kinsman's seas, the variations are slight as 
is shown by the following comparison: 

Kinsman - record 072 OAWF - 3 Fans 

Skewness +.092 -.02*15 

Kurtosis -.031 -.35*»9 

The mean is also seen to be non-zero but is less 
than 0.1 centimeter. 

3. Autocorrelation Function 

The autocorrelation function of the ocean wave 
heights is characteristically that of band limited noise and 
takes the shape of an attenuated cosine [Ref. 8]. In finding 
the correlation function of the OAWF three fan ocean, 
essentially the same setup was used as in the previous two 
sections. However the output from the wave probe, 
characterizing the surface, was terminated, not in the 
OPHELEA system as before, but in a PAR CORRELATOR. This 
device is capable of computing either auto- or cross- 
correlation functions and displaying the result on an 
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Figure 17. W&ve Height Histogram and Corresponding Normalized 
Gaussian PDF 



oscilloscope or a recorder. In this case an X-Y RECORDER 
was used. A diagram of the experimental setup is shown in 
Figure 18. 

An average surface wave height correlation function 
for time lags, x = 0.0 to 2.0 seconds, was arrived at by 
superimposing the eleven runs printed on the X-Y RECORDER 
and shown in Figure 19. The resultant average is shown in 
Figure 20 and can be seen to follow the expected form. The 
rate of "attenuation" is initially somewhat les^ than 
measured by Kinsman, with the OAWF correlation at the first 
peak down to 0.7 compared to 0.50 for Kinsman's ocean data. 
The second peak correlation conversely has dropped to about 
0.3 while Kinsman's C(x) remains at about 0.4 after time 
lags of two ocean periods. 

The correlation function loses its consistency after 
two ocean periods as may be seen from Figures 19 and 21. 
Figure 21 shows the function computed with a precomputation 
period = T = 2.0 seconds set into the PAR CORRELATOR. With 
this setting the correlation function is computed for 
x=2.0 to 4.0 seconds or in effect the function of Figure 19 
is continued for another two seconds. For this longer time 
lag, it is difficult to see any repeatability in successive 
correlation functions. 

While the OAWF values are somewhat different from those 
of the Kinsman ocean, the similarity in shape is obvious. 
This validation of the ocean model will be of importance 



41 




Figure 18 . Experimental Setup for Determining the OAWF 
Surface Wave Autocorrelation Function 
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Figure 20. Average OAWF Surface Wave Height Autocorrelation Function 
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later in this work when comparisons are made with the 
modulation envelope of scattered sound. 



4 . The Roughness Parameter 

The previously mentioned roughness parameter , -\fg , 
is defined by Beckmann and Spizzichino [Ref. 2] as: 



yg = [cos 0^ + cos e 2 ] 



with 



a 

X 

0 1 J0 2 



ocean RMS wave height 
acoustic wave length 

angles of incidence and reflection, respectively 
(measured from the normal to the ocean surface). 



It is easily seen that a ’’smooth" condition (VF << 1) can 
be achieved through use of low acoustic frequencies or large 
angles of incidence (shallow grazing angles) as well as with 
the obvious small RMS surface wave heights. The converse is 
of course also true for "rough" conditions >> 1). 

VF was the factor that was used to model the OAWF ocean to 
the actual ocean in later research and an example of typical 
sets of conditions in both domains yielding the same 
-\Jg = .076 is given below: 



Wind 

RMS wave height 
Acoustic Frequency 
Source/Rcvr Depth 
Source/Rcvr Separation 



OAWF 


OCEAN 


3 Fans 


20 KTS 


. 7 cm 


54 cm 


13 kHz 


300 Hz 


10 cm 


200 m 


150 cm 


7000 m 



As noted above i/F is the same for both sets of 
parameters and thus the typical ocean conditions noted in 
the right column can be reproduced if the values in the left 
column are employed in the OAWF. 

B. COMPARISON OF WAVE HEIGHT AND DEMODULATED SCATTERED 

SOUND SPECTRA AND CORRELATION FUNCTIONS 

As noted in the Theory section the frequency spectrum of 
the amplitude modulation of the scattered sound is directly 
related to the frequency spectrum of the wind driven waves. 
Hence it is important to compare these frequency spectra. 
Further it is of interest to compare autocorrelation functions 
in order to gain an insight into signal processing possibili- 
ties. In addition to these two comparisons the cross 
correlation function was computed. 

1. Frequency Spectra Comparison 

Though some trials described later in this work were 
conducted with a directional source, the principal study 
used an omnidirectional source and receiver, which were felt 
to best model the submarine/sonobuoy scenario. The source 
was a locally-constructed transducer made from a 2" (o.d.) 
Glennite Ceramic sphere which was covered with several coats 
of neopreme . This source transducer will be referred to 
hereafter as "OMNI". The receiver was a standard Atlantic 
Research LC-10 hydrophone and its response was immediately 
amplified by 30 dB, NUS Corporation FET preamp powered by 
a 12 volt DC POWER SUPPLY. Both the LC-10/NUS combination 
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and the OMNI, as well as a directional piston type transducer 
(P-27) used later, are shown in Figure 22. 

As the spectra of the OAWP sea waves had already 

I 

been computed from an average of forty runs, only the 
average sound amplitude modulation spectrum remained to be 
determined. To accomplish this the scattered AM sound was 
demodulated by the detector shown in Figure 14 and the 
resulting time-varying DC signal sent to the OPHELEA system. 

A photograph of the wave forms of the AM scattered sound, 
the demodulated scattered sound and the instantaneous wave 
height are shown in Figures 23 and 24. A complete schematic 
of the experimental setup is shown in Figure 25. 

Forty runs were conducted with 256 A/D samples of 
the demodulated acoustic 20 kHz signal being taken for each 
run at a sampling frequency of 20 Hz, again using program 
PING 1C. As before, a 128-point spectral density versus 
frequency readout was calculated for each run and the 
resulting forty spectra averaged and normalized for comparisons 
with the previously determined average wave spectrum. These 
two averaged, normalized spectra are plotted in Figure 26 
which shows the obvious similarity. The demodulated sound 
spectrum has a shape that is remarkably similar to its 
counterpart though its peaks are about 0.2 Hz less than the 
predominant wave frequencies. This difference is presumably 
due to the surface drift velocity [Ref. 9]- The 3 dB bandwidth 
of the demodulated sound spectrum is about 0.5 Hz as compared 
with 0.6 Hz for the wave spectrum. 
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Figure 22. Hydrophone (LC— 10) and Transducers Used 
in this Research 
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Figure 23. 

Amplitude Modulated 
Sound (top) and 
Instantaneous Wave 
Height 



Figure 2b . 

Demodulated AM Sound (top) 
and Instantaneous Wave 
Height 
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Figure 25. Experimental Setup for Determining the Spectrum 
of Demodulated Scattered Sound 
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;ure 26. Normalized Frequency Spectra of Demodulated Forward Scattered Sound 
and OAWF Surface Wave Height 



These frequency spectra, though exhibiting some minor 
differences, are considere d"'t o - b j H.-trs sentially the same and ' 
confirm the theoretically expected result. 

' 2 . Autocorrelation Function Comparison 

The same procedure was used in determining the 
demodulated scattered sound autocorrelation function as in 
the case of the wave height version. The output of the 
previous setup was terminated in the PAR CORRELATOR instead 
of OPHELEA and recording was done on the X-Y RECORDER as 
before. Again eleven runs were conducted for delays 
x = 0.0 to 2.0 seconds. These recordings are shown in 
Figure 27. An average autocorrelation function was then 
fitted to this record. A plot of both records, wave height 
and demodulated scattered sound, is shown in Figure 28. 

Though the ordinate at t = 0 is slightly different for the 
two functions due to idiosyncracies of the PAR CORRELATOR 
and X-Y RECORDER, it is apparent that the results are nearly 
identical for the two cases. Further, it is evident that 
significant correlation for both variables exists for three 
"peaks", that is, over two ocean periods. The demodulated 
sound correlation function also becomes confused for delays 
greater than 2.0 seconds, as is the case with the wave height, 
as seen in Figure 29. It should be noted that the vertical 
scale in Figure 29 is expanded over that of Figure 21 and no 
inference of greater correlation should be drawn. 
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Figure 28. Average Autocorrelation Functions of Demodulated 
Forward Scattered Sound and Wave Height 
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3. Cross Correlation Function 



Using the two setups j. be d in the autocorrelation 

section above the two wave forms were cross correlated, again 
using the PAR CORRELATOR and X-Y RECORDER. The same sort of 
averaging process was then performed on the record shown in 
Figure 30 and yielded an average cross correlation function 
which is shown in Figure 31. The t = 0 ordinate varied 
considerably between positive and negative values of 
correlation and the choice of a near zero initial value was 
purely arbitrary. Again three peaks are seen to dominate, 
at about t = 0.8, 1.25 and 1.7 seconds. This is to be 
expected from the cross correlation function of two signals 
which both exhibit significant correlation over two ocean 
cycles (or have three "peaks" in their autocorrelation 
functions ) . 

C. EFFECT OF FREQUENCY AND ANGLE OF INCIDENCE 

ON SCATTERED SOUND MODULATION 

To determine whether the modulation frequency of the 
scattered sound remained essentially constant with changing 
acoustic frequency and varying angle of incidence, a series 
of trials was conducted with different combinations of these 
two parameters. It was felt at this stage of the research 
that there might be some advantage in using a directional 
hydrophone to make the detection. 
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1. Changing Frequency and Angle of Incidence 

In these tests a directional transducer, a USRD 
Type F-27, was mounted in a trainable bracket which permitted 
angles of incidence, measured from the surface normal, of 
from 0° to 90°. The F-27 is shown in Figure 22. The 
scattered sound was received by an omnidirectional LC-10 
hydrophone which was positioned at the same depth as the 
F-27 and 2.75 meters away. To achieve different angles of 
incidence while retaining specular scatter the F-27 and 
LC-10 were positioned at depths of 0.99, 0.^9 and 0.27 meters 
which equate to angles of incidence of approximately 58°, 

70° and 79° respectively. Prior to each trial the F-27 was 
elevated and depressed, with the surface quiescent, until 
the maximum response ensured that the observation was of 
specular scatter. Ten runs, using program PING 1C (Appendix 
A), were conducted at each position with 256 A/D samples 
being taken at a sampling rate of 20 Hz. A 128-point spectral 
density versus frequency output was generated as in the 
determination of previous spectra. Again the WAVETEK 
oscillator provided the pulse train for proper A/D converter 
triggering at low frequencies. A diagram of this setup is 
shown in Figure 32 and a table of run parameters and 
particular equipment settings is found in Table II. 

After the above 80 runs were completed, the ten 
spectra from each angle/frequency combination were averaged 
and normalized to permit comparison. The results were then 
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Figure 32. Experimental Setup for Determining the Effects 
of Varying Angle of Incidence and Frequency 
on Scattered Sound Modulation 



ill 



plotted as Indicated in Table I below with results compared 
for constant angle of incidence and variable frequency and 
vice versa. 



TABLE I 

•SUMMARY OF VARIABLE ANGLE OF 
INCIDENCE AND FREQUENCY TRIALS 

Figure Acoustic Frequency Angle of Incidence 



33 


30,40 kHz 


79° 


3 |4 


20,30,40 kHz 


O 

O 

c*- 


35 


20,30,40 kHz 


0 

CO 

in 


3 6 


20 kHz 


70° ,58° 


37 


30 kHz 


79°, 70°, 58 


38 


40 kHz 


79° , 70° ,58 



As may be seen from Figures 33 through 38 the 
frequency spectral peaks differ by less than 0.3 Hz in all 
cases except Figure 35. The somewhat wider spread, about 
0.5 Hz, between peaks is thought to be due to greater 
contributions from the side lobes at that angle of incidence. 
Even this difference, though larger, is comparable to the 
3 dB bandwidth of the wave frequency spectrum. 

Based on this work, it was concluded that the 
modulation spectra of the scattered sound received by a 
point hydrophone are similar for variations in acoustic 
frequency or angle of incidence. It is felt that the 
differences are due to inadequate averaging and/or the 
dependence on side lobes. 
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Summary of Parameters Affecting Trials with Variation 
of Frequency and Angle of Incidence 
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Figure 33. Normalized Frequency Spectra of 30 kHz and 40 kHz Demodulated 
Sound at an Angle of Incidence of 79° 
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Figure 3’i. Normalized Frequency Spectra of 20 kHz, 30 kHz and 40 kHz 
Demodulated Scattered Sound at an Angle of Incidence of 70 
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Figure 35- Normalized Frequency Spectra of 20 kHz, 30 kHz and 40 kHz 
Demodulated Scattered Sound at an Angle of Incidence of 58 
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Figure 36 . Normalized Frequency Spectra of 20 kHz Demodulated Scattered Sound 
at Angles of Incidence of 38° and 70° 
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Figure 37. Normalized Frequency Spectra of 30 kHz Demodulated Scattered Sound 
at Angles of Incidence of 58°, 70° and 79° 
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Figure 38 . ’formalized Frequency Spectra of 40 kHz Demodulated Scattered Sound 
at Angles of Incidence of 58 ° , 70° and 79° 



2 . Modulation Strength with Changing Receiver Angle 

The object of this set of trials was to determine if 
the modulation effect could be increased by training a 
directional receiver to a preferred angle of reflection with 
a point source. Although acoustic reciprocity would be 
expected to hold, in this case control would be at the option 
of the listener. Specifically it was felt that possibly if 
the receiver hydrophone were directed normal to the modulating 
ocean surface, the relative modulation of the received 
scattered acoustic signal would be increased. The experi- 
mental setup was essentially as before; however, now the 
P-27 was the receiver and the previously referred to OMNI was 
the source, with both at a depth of 1.0 meters. The three 
fan OAWP sea was again the environment with the water level 
standardized as before. A setup diagram is shown in Figure 
39. Ten runs were conducted at each of seven angles of 
incidence ranging from direct path, through specular to 
normal incidence. The acoustic frequency was 30 kHz and 
amplifier and filter settings were kept constant for all 
runs so that valid comparisons could be made. Table III 
shows a summary of these trials. 

Each of the seventy runs in Table III consisted of 
taking 256 A/D samples at a 20 Hz sampling frequency (using 
program PINC 1C) as before; then averaging the spectral 
density information for each angular setting. The plot of 
these average spectra is presented in Figure *10 which shows 
log spectral density plotted versus scattered sound 
modulation frequency. 
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TABLE III 



SUMMARY OF VARIABLE -45K?;,F- .RBCEJYSR TRIALS 



Run Numbers 



Angle of Incidence 



61-70 

71-80 

81-90 

91-100 



45 ° 

30° 



90° (Direct) 
58° (Specular) 



101-110 

111-120 

121-130 



0° (Normal 



25 ° 

75 ° 



Note: Again the angle of incidence is measured 

normal to the surface. 



All levels have been plotted without adjustment 
so that the amplitude of the modulation at different angles 
of incidence can be seen simultaneously. Summarizing: 

a. The greatest modulation strength was observed 
at 0 = 58°, that is the specular direction. 

b. Modulation amplitude declined regularly as the 
angle of receiver is changed in either direction from the 
specular. 

c. Although the incoherent component (modulated 
amplitude) may possibly be increased relative to the 
coherent component, the modulation amplitude is down 30 dB 
(re specular) with the receiver trained normal to the sea 
surface (0 = 0° ) . 

Knowing that the amplitude modulation is greatest at 
or near the specular direction might provide a means of 
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arriving at a rough fix or target probability area. Assuming 
that target bearing could be provided, by a directional hydro- 
phone and that the local sound velocity profile is known, the 
approximate range could be found by applying ray path 
techniques once the angle of greatest modulation (specular 
direction) is found. This scheme would of course require 
a hydrophone system which is directional in both azimuth 
and elevation. 
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V. ENHANCING THE TIME-VARYING SIGNAL 



At this stage, with reasonable assurance that the OAWF 
sea accurately modeled the actual ocean and that the modula- 
tion spectrum of the scattered sound was essentially that of 
the surface waves, efforts were begun toward improving signal 
processing using the model and the information gained to date. 
If the ocean wave spectrum were of a single fixed frequency 
the sound would fluctuate at that same single fixed frequency 
and sound signal enhancement would then be easy. It would be 
necessary merely to identify a loud point in the received 
signal and then to listen at later intervals equal to the 
ocean wave period. The ocean is regrettably not this simple 
as is obvious from the respective frequency spectra and 
therefore alternative or modified processing means had to be 
identified and implemented. 

A. PROGRAMS PERK 1A AND PERK IB 

To begin to explore the problem, a pair of computer 
programs modestly named for this. rookie acoustician were 
developed which would accomplish the following: 

1. PERK 1A orders the computer to direct the A/D 
converter to take a group of samples (say 6^), hereafter 
called a "BLOCK", and then to "pause". 

2. The "pause" is triggered by the computer taking its 
last ordered sample and PERK 1A then shifts the computer 
into a Fortran "Do Loop" which has an index proportional) to 
the desired time delay. 
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3. When the "Do Loop" has completed Its preset number 
of operations (and the delay or pause has been completed) the 
program signals the A/D to take another BLOCK of samples. 

This Is followed by another "pause", another BLOCK, etc. 

The net result is a series of BLOCKS of A/D samples 
taken at specified intervals along the modulated signal 
which, when a total of 8192 data points have been taken (the 
present storage limit of OPHELEA), are read out to a tape 
cassette. This data tape is then put aside as the second 
program, further modestly entitled PERK IB, is loaded. 

PERK IB reads the time series data of the previously 
recorded sample BLOCKS back into memory destructively and 
performs an FFT spectral analysis on each BLOCK. The program 
may then be directed to output all of the complete spectral 
analyses, or just the major spectral peak in each analysis. 
Thus a sequence of successive sound spectral density peaks 
is available, each separated by a fixed time delay, the 
"pause". The duration of the "pause" is selected to provide 
the desired number of spectral densities of predetermined 
frequency resolution during an "ocean period". An "ocean 
period" is here defined as the reciprocal of the peak 
frequency in the ocean spectrum. The effect of the modulation 
caused by the ocean is vividly illustrated by Figure ^1, which 
depicts the peak sound spectral density versus time, and which 
shows a swing of about 18 dB between maxima and minima of 
signal strength. Copies of PERK 1A and PERK IB may be found 
in Appendix A. 
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Figure 4l. Changing Spectral Density Amplitude of a Forward Scattered Acoustic 
Signal. (Weak Signal - Large Fluctuations) 



In this work, the BLOCKS almost invariably contained 
64 A/D samples and the "pause" was selected to provide eight 
BLOCKS per ocean period. The 'following is a sample delay 
time calculation for a 20 kHz acoustic signal, modulated by 
a 2.5 Hz ocean (ocean period 0.4 sec.) and sampled at 80 kHz. 

For 64 Samples/BLOCK at 80 x 10^ samples/sec.: 
sampling time/BLOCK =0.8 msec. 

For 8 BLOCKS/ocean period: sampling time/ocean period = 

6.4 msec. Therefore, for one ocean period, total delay time 
= 400 msec - 6.4 msec = 393*6 msec. 

Delay between BLOCKS = 393.6/8 = 49.2 msec. 

Total data collected in 16 ocean periods 

_ / 64 samples >, , 8 BLOCKS , / l6 ocean periods x 

^ BLOCK ocean period^ run ' 

_ g -,02 samples 
y run 

For the purposes of this research eight BLOCKS per ocean 
period were considered adequate coverage. Each of the eight 
BLOCKS was said to be at a different PHASE lettered A through 
H. Every eighth BLOCK would be at the same phase. Again if 
the ocean were monochromatic and the spectral density peaks 
in each PHASE averaged, one average, corresponding to the 
strongest part of the signal, would predominate and the other 
PHASES of lesser amplitude would be distributed sinusoidally 
around it . 

PERK 1A and PERK IB are highly versatile programs with a 
number of options available as summarized in the following 
table (Table IV). 
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TABLE IV 



SUMMARY OF PERK 1A AND PERK IB PARAMETER OPTIONS 



Number of data points/BLOCK 



Total number of A/D data points 



up to 8192 

from 16 to 512 in powers 



of 2 (2 N ) 



Delay between BLOCKS 



0.01 to 999.9 msec with 
accuracy to 0.01 msec. 

2 to 16 in powers of 2 



Number of BLOCKS/ocean period 



= number of PHASES 

* 

The program also provides a unique number for each run, 
based on day, month and run sequence of that day (i.e., 

#15101 is the first run conducted on 15 October) , and also 
prints out input parameters such as sampling frequency, prior 
to displaying the output. 

It should be noted at this point that PERK 1A/1B 
were tested on a known sinusoid of 2.5 Hz and found 
satisfactory. In that test, groups of 8 A/D samples separated 
by 0 . 4 seconds were taken at 80 kHz and plotted versus time. 

No visible discrepancy existed over 24 sinusoidal periods. 

To recapitulate, the programs PERK 1A and PERK IB were 
used to take BLOCKS of 64 A/D samples, properly spaced 
temporally to provide eight PHASES (A-H) per ocean period, 
over sixteen ocean periods. This equates to 128 BLOCKS or 
16 BLOCKS per PHASE. Samples were taken at four times the 
acoustic frequency, or twice the Nyquist rate, to diminish 
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the effects of aliasing. A pictorial representation of 
this signal processing scheme is shown in Figure 42. 

B. EFFECT OF SIGNAL STRENGTH ON MODULATION - 

LLOYD MIRROR MODELING 

One of the most obvious times when signal processing 
enhancement is important is when the received acoustic 
signal is weak and may be fluctuating in and out of the 
background noise. In order to provide a common basis for 
comparing varying signal strengths, it was decided to employ 
different aspects of the Lloyd mirror phenomenon to observe 
varying strengths of the output signal. 

1 . Brief Review of Lloyd Mirror Effect 

As every school-boy acoustician knows, the Lloyd 
Mirror Effect is caused by constructive and destructive 
interference between the direct and smooth surface reflected 
paths taken by an acoustic signal. For a mirror surface within 
a Critical Range R c = , the effect is one of pressure 

peaks and valleys, with the peaks falling off as 1/r and the 
valleys going to zero (h = receiver depth, d = source depth, 

X = acoustic wavelength). Beyond the critical range the 

p 

pressure declines as 1/r . The OAWF tank is capable of 
measurements both within and beyond the critical range for 
frequencies in the vicinity of 13 kHz. Another means of 
visualizing this effect is by replacing the reflected rays 
by an image source at I as also shown in Figure 43. 

In the case of a rough surface the peaks caused by 
reinforcement arc not as great and the coherent scattering 
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Figure >\2. PERK 1A/PERK IB Signal Processing Scheme 
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Figure ^3. Lloyd Mirror Effect 




Figure M . Differences in Lloyd Mirror Effect for 
Rough and Smooth Surfaces 
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term, given by the first Dirac delta term in the expression 
for the scattered sound spectrum (Theory section), predominates. 
Further the valleys caused by cancellation are now minima 
instead of zero as the incoherent or modulated term in the 
S(J2) expression becomes important. This effect is shown in 
Figure 44 which contrasts the rough and smooth situation. 

The setup used in the following three trials was essen- 
tially the same as in all signal processing runs. The OMNI 
was the source and the LC-10 served as receiver (both at a 
depth of 10 cm.), with a 3 Fan OAWF sea, set to the standard 
water level, modulating the transmitted 20 kHz acoustic 
signal. A/D samples were taken with program PERK 1A at a 
sampling rate of 80.0 kHz and analysis was performed by 
PERK IB. A setup diagram is shown in Figure 45. The rack 
on which the LC-10 was mounted was moved back and forth with 
the water surface quiescent (fans off) until a Lloyd mirror 
maximum and adjacent minimum were identified. A third 
position was then chosen arbitrarily between the two extremes 
and designated "Intermediate". Runs were conducted at each 
of these three locations and a summary is presented in Table V 
below. 

From the output of PERK IB it was possible to compare the 
signal fluctuations at the three positions. At the Lloyd 
mirror maximum, corresponding to a strong signal (larger 
constant amplitude coherent and smaller variable incoherent 
contribution from the surface scattered sound), the maximum 
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TABLE V 



TRIALS CONDUCTED AT VAET0U5,.,J LCJfH MIRROR POSITIONS 



Run No. 


Lloyd 

Mirror 


Samples/ 

BLOCK 


Delay No. of 

PHASES 


Total Ocean 
Periods Covered 


10101 


Min 


64 


47.69 8 

msec 


16 


10102 


Max 


64 


(assumed 8 

ocean 
freq. 


16 


10.03 


Inter- 

mediate 


64 


2.578 Hz) 8 


16 



spectral density peak exceeded the overall average by about 
2.1 dB and the minimum peak fell about 2.0 below the average. 
The total range of amplitudes was about four decibels. 
Applying the same comparison technique to the other two runs 
yields a swing of about 8.5 dB at the "Intermediate" position 
and almost 19 dB at the position of the weak signal where 
there is a maximum fluctuation of the amplitude. 

It is apparent that in this model fluctuations do exist 
regardless of signal strength and that they are of greatest 
importance when the signal is weak. This further points up 
the importance of adapting knowledge of the sound modulation 
to enhancing a signal processing method as a means of 
acquiring and/or retaining low strength signals in a noisy 
or highly variable environment. 





80 kHz 



-> 



TO OPHELEA 
(Sampling Freq.) 



Figure A 5. Experimental Setup for Determining the Effect 
on Spectral Density Variation with Time at 
Three Lloyd Mirror 1 Positrons 
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C. "BEST PHASE" METHOD 

As previously mentioned, ,±f t-.he.. scattered sound were 
modulated at a single fixed frequency there would be a PHASE 
which, when once correctly chosen, would always yield the 
strongest average signal. The modulation frequency is of 
course not monochromatic as shown earlier in Figure 4l and 
in fact its spectrum has about a 0.5 Hz half-power bandwidth 
(Figure 26). 

This not withstanding it was felt both instructive and 
possibly rewarding to compare the relative strengths of the 
eight PHASES. Thus an annex was appended to PERK IB which 
averaged all of the spectral density peaks in each PHASE 
(every eighth peak in this case) and printed out a PHASE 
average as well as the overall average of all peaks. These 
values for 8 runs are shown in Figures 46(a-h), which depict 
PHASE average spectral density plotted versus PHASE with the 
dashed line indicating the overall average. All of the runs 
are seen to, in fact, follow a rough sinusoidal pattern 
particularly Figures 46(c) and 46(g) but others such as 
Figure 46(a) appear to be at double the ocean frequency. 

The Figures 46(a~h) show spectral density on a linear scale. 

If the "BEST PHASE", or highest PHASE average is compared 
to the overall average, the gain is in the range of 0.3 dB 
(for the Lloyd mirror max-strong signal case) to 0.5 dB 
(for the weak signal), which is hardly anything to write 
home about. 
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AVERAGE SPECTRAL DENSITY (Linear) AVERAGE SPECTRAL DENSITY (Linear) 




(a) Run 18901 




PHASE 

(b) Run 20901 

Figure J l6. Average PHASE Spectral Densities 
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AVERAGE SPECTRAL DENSITY (Linear) AVERAGE SPECTRAL DENSITY (Linear) 
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PHASE 
(c) Run 20902 




PHASE 
(d) Run 10101 

Figure *l6. Average PHASE Spectral Densities 
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AVERAGE SPECTRAL DENSITY (Linear) AVERAGE SPECTRAL DENSITY (Linear) 




Overal 1 
Average 



C D E 
PHASE 

(e) Run 10102 



26 

25 

24 

23 

22 

21 

20 



1 



Overal 1 



T 



Average 



A B C D E F G II 
PHASE 

(f) Run 10103 

Figure 46. Average PHASE Spectral Densitie: 
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AVERAGE SPECTRAL DENSITY (Linear) AVERAGE SPECTRAL DENSITY (Linear) 




PHASE 
(g) Run 15101 




PHASE 
(h) Run 15102 

Figure ^6. Average PHASE Spectral Densities 
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In summary the "BEST PHASE" method applied innocently 
over 16 ocean periods is seen to offer some small advantage. 
However, because of the broad band modulation of the scattered 
sound modulation and the drift of the best phase over many 
ocean periods the gain realized is minimal. 

D. "MONDAY MORNING SIGNAL PROCESSING" METHOD 

The next approach began with hand plotting the spectral 
density versus time information presented by PERK IB and 
examining it to see if any combinations would become apparent. 
Samples of the plots are included in Figures 4l, 47 and 48 
and show the fluctuations of spectral density with time for 
the Lloyd mirror minimum, "intermediate" and maximum, 
respectively. 

To analyze this displayed information a clear plastic 
mask was constructed with a sinusoid of the same frequency 
as the assumed peak OAWF ocean frequency printed on it. The 
mask was then moved along the time axis of the spectral 
density versus time plot to determine how many peaks or near 
peaks corresponded to the ideal sinusoid peaks. Typically 
groups of at least three but up to five woul d line up prior 
to an irregularity or "glitch" which would require shifting 
phase to a subsequent group. This was not unexpected as the 
auto- and cross-correlation functions discussed earlier in 
this work indicated good correlation for at least two or 
three ocean periods. 
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Figure ^7- Changing Spectral Density Amplitude of a Forward Scattered Acoustic 
Signal ("Intermediate" Strength Signal) 
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Figure 3. Changing Spectral Density Amplitude of a Forward Scattered Acoustic 
Signal (Strong Signal - Small Fluctuations) 



The above procedure was felt to simulate a signal 
processor which would Identify a strong signal and then 
sample later, only at Intervals which are multiples of the 
ocean period, three to five times to gain signal enhancement 
over continuous Integration. The results of this signal 
processing In hindsight — hence the Monday Morning label — 
are Interesting and are Illustrated below: 

Referring to Figure 4l, the Lloyd mirror weak signal 
(large fluctuations) case; If the BLOCKS indicated by arrows 
are chosen and averaged as follows : 



Group 

1 

2 

3 



BLOCKS 

5,13,21 

39,^7,55,63, 71 
98,106,114,122 



M 2 SP Ave. 
‘ '44.70 
25.20 
42.55 



These groups of BLOCKS correspond to identifying initial 
strong signal (i.e. BLOCKS 5,39 or 98), then sampling at the 
ocean frequency for only three to five peaks further. The 

payoff becomes apparent when these Monday Morning averages 
are compared to the overall average of 19-6 which yields 
gains of +3-6, +1.1 and +3-4 dB for groups 1, 2, and 3, 
respectively . 

Gains in the case of the modeling strong signal (weak 
fluctuations) ranged from +1.0 to +1.4 dB and those for the 
Lloyd Mirror "Intermediate" response were in the +0.4 to 
+1.1 dB range. The interested reader is invited to verify 



these values, or try to improve on them with other 
combinations, by marking a ruler or like object in intervals 
of 0.4 seconds and sliding it along Figures 4l, 47 and 48. 

. The advantage of Monday Morning Signal Processing, is 
obvious in the weak signal environment (where it is most 
needed, of course) and is at least worthwhile, though not 
vital, when observing a stronger signal. The author will 
consider modifications to this technique and suggestions for 
real time analysis in the final section. 



E. SUCCESSIVE AVERAGING 

In an attempt to validate statistically the "seat of 
the pants" notion that three to five ocean periods is the 
maximum extent to which sampling at multiples of the ocean 
period can be profitably employed, a system of averaging was 
developed and appended to PERK IB (which became PERK 1C). 

The spectral density peaks were first divided into eight 
PHASES (A through H) of sixteen peaks each. These PHASE 
groups were then blindly averaged in ensemble group as follows: 

Taking PHASE B for example. 



Average by two's 



B1+B2 B2+B3 

2 3 3 



B16+B1 

o 



Average by three's 



B1+B2+B3 B2+B3+B4 

3 3 3 



B16+B1+B2 

3 



et c . 

Up to Average by eights 
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These averages were then compared with the overall 
average and the number that exceeded it noted. Rather than 
confirming the initial hypothesis, however, 'the results were 
contradictory; indicating that groups of three to five were 
most profitable in some PHASES while averages of five to 
eight were most beneficial in others. In fact one PHASE 
indicated that it made little difference how many were 
averaged — that all combinations were equally profitable. 

The author remains convinced however, that, though it 
has not been statistically proven, when compared with blind 
averaging, the Monday Morning technique demonstrates that 
averaging more than three to five peaks in a group is 
counterproductive . 

In summary the blind Best Phase method provides only a 
minimal gain of about half a decibel when applied to a weak 
signal while the Mondary Morning technique provides nearly 
a +3 dB enhancement of the signal processing for very weak 
signals if sampling is limited to three to five ocean cycles. 
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VI . FUTURE RESEARCH, CRYSTAL BAL L GAZING AND OTHER MUSINGS 

The author feels that the research described in the 
foregoing pages has only scratched the surface of a viable 
means of acoustic signal enhancement. Some basic parameters 
and features of the modulation of surface scattered sound 
have been identified or confirmed and in some cases loosely 
quantized. These ideas may be of immediate use to signal 
processors, however additional research is felt to be 
indicated in the following two areas: 

A. Further investigation of the phenomenon for additional 
clues which might increase its value/utility and prove more 
rigorously some notions presented herein. 

The initial concept of identifying a signal peak (or loud 
point), then sampling later at multiples of the ocean period, 
which at this point might be charitably described as 
simplistic, could be more effective in a modified form. If 
a peak were identified and the next sampling interval were 
sought within the interval T±AT, where T is the ocean period 
and AT is related in some unknown way to the bandwidth of 
the ocean spectrum and/or the rate of "attenuation" of its 
autocorrelation function, greater improvement might be 
realized. Of course the interval of sampling would increase 
by additional multiples of AT if the starting point were not 
reset at a peak after each group of samples was taken and 
the peak identified. (See Figure J )9.) It is felt that the 
previously mentioned limit of three to five sampling intervals 
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Interval not reset to signal peak; 
next sampling interval must be 
±2AT to cover all possibilities. 
Then ±3AT, etc. until continuous 
integration results 
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Figure ^9. Improving on the "Best Phase" Processing Method 



before the diminishing returns point is reached may also 
apply here . 

It might be even more profitable to take samples at a 
rate of several BLOCKS per ocean period over some observation 
interval (greater than the ocean period) as done in this 
research and then to compute, sort and rank the spectral 
density peaks by magnitude. Once so arranged the following 
options would be available: 

1. use only the highest peak in each observation interval 

2. use only the highest peak in each ocean period and 
average over the observation interval 

3. 'use the five (or any other arbitrary number) highest 
peaks observed and average as in 2.). 

4. use only those peaks which exceed the overall average, 
which would also have to be computed, by 3 dB ( or 
employ some lesser standard, such as 2 or 3, if none 
did) . 

As an aside it should be noted that when speaking of the 
"highest spectral peak" it is implicit that the highest for 
each frequency observed is computed, compared and utilized. 
This research was conducted with a single acoustic frequency 
but future trials might include several tranducers operating 
at various frequencies and signal output levels. A BLOCK of 
more than 6 l 1 samples may be required in this case to provide 
the desired frequency resolution. 

Another possibility might be to continuously process for 
most of the "-loud" half of the ocean period (5 seconds or 
so in the actual ocean) , then process and use the received 
frequency spectral information to image correlate during 
the quieter half cycle. Again the three to five period 
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profitable limit is felt to pertain. This method might 
provide even better results obtained by active 

sonar image correlating as the changing doppler effects and 
attendant correlation mismatches would not be a factor over 
such short time spans. 

B. The second major area of suggested study would be to 
spectrally analyze, plot, and employ the Monday Morning Signal 
Processing Technique on actual passively-recorded submarine 
tapes to see if gains comparable to this research are realized. 
Additionally the modulation frequency spectrum and autocorre- 
lation functions could be determined and compared with both 
the locally observed wind/sea data and the results of this 
work. 

In summary, there is no doubt as to the fact that sound 
scattered from the ocean surface is amplitude modulated at 
the basic surface wave frequencies. The ability to estimate 
the ocean wave frequency and spectrum is of course vital to 
taking advantage of this phenomenon and has not been addressed 
to this point. The author feels that ballpark estimates 
could be arrived at on station through a knowledge of fetch 
and surface wind speed and duration and assuming a Pierson- 
Moskowitz wave height spectrum. Further it is felt that 
meteorological organizations such as Fleet Numerical Weather 
Central could supply this information if tasked to do so, 
so that signal processors would be provided with up-to-date, 
tactically-useful data. 
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Though other mechanisms such as internal waves and tem- 
perature microstructure also contribute to fluctuations, the 
surface scattering is of particular importance for near- 
surface sources and/or receivers. What remains in the author’s 
opinion is to find the means to most effectively use the 
various ocean modulation effects to enhance signal processing. 
It is felt that clever signal processors, if properly in- 
formed, could devise numerous techniques far more subtle and 
comprehensive than the rather simplistic methods suggested 
in this and the foregoing sections and that even greater 
gains might be realized. It should be emphasized that use 
of these techniques is not contended to be the miracle cure 
for the Anti-Submarine Warfare problem. The author does feel 
however that ocean modulation phenomena, if used to best 
advantage, will provide part of the answer required to cope 
with quieter next generation targets. 
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